Spin relaxation and combined resonance in two-dimensional 
electron systems with spin-orbit disorder 
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Disorder in spin-orbit (SO) coupling is an important feature of real low-dimensional electron 
structures. We study spin relaxation due to such a disorder as well as resulting abilities of spin 
manipulation. The spin relaxation reveals quantum effects when the spatial scale of the randomness 
is smaller than the electron wavelength. Due to the disorder in SO coupling, a time-dependent 
external electric field generates a spatially random spin-dependent perturbation. The resulting 
electric dipole spin resonance in a two-dimensional electron gas leads to spin injection in a frequency 
range of the order of the Fermi energy. These effects can be important for possible applications in 
spintronics. 

PACS numbers: 72.25.Rb, 72.25.Hg 
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Electron dynamics in low-dimcnsional semiconductor 
structures reveals features of a spin dependent trans- 
port that are interesting for fundamental and applied 
researcbl. One of the main ingredients necessary to gen- 
erate spin dependent transport in nonmagnetic semicon- 
ductor systems is the SO interaction. Such an inter- 
action offers a possibility of an efficient and fast spin 
manipulation with electric fields, which in turn allows 
to prepare a required spin statei^i^i^i^iii^ At the same 
time, spin relaxation and decoherence due to the SO cou- 
pling prevent long-distance spin propagation. Two mod- 
els are widely used to describe the SO coupling in low- 
dimensional structures: the Rashba and the Dresselhaus 
ones. In both models, the SO field and the correspond- 
ing spin precession rate are approximately linear in the 
electron momentum. Random evolution in the momen- 
tum due to collisions with impurities, phonons, and other 
electrons results in randomness in the spin precession, 
and thus leads to spin relaxation. However, in reality 
both interactions have an intrinsic randomness due to 
system imperfections, including the fluctuations in the 
dopant ion density2.iiS or random bonds at the quantum 
well (QW) interface^^. Even if the mean values of the 
Rashba and Dresselhaus fields vanish, their fluctuations 
remain and can cause interesting consequences, includ- 
ing memory effects^, spin Hall effect in the finite-size 
systemsi^.^ and spin-dependent localization.''^'^ 

There are at least four different two-dimensional (2D) 
systems, where the SO disorder plays an important or 
crucial role. First, the effect of random SO coupling can 
be responsible for the spin relaxation in Si/Ge QWs.^^'JA 
Second, the spin-dependent disorder influences^ the 
spin helix pattern recently observed in the GaAs (001) 



QW with the balanced Rashba and Dresselhaus termSf^^ 
Third, the randomness causes relaxation of the s pin 
component along the growth axis observed in Ref.[l7[ 
in GaAs (Oil) QW, investigated now for spintronics 
applications. ^^'^'^ Fourth, the most recent example of the 
system with random SO coupling is graphene, where the 
randomness and spin relaxation appear due to the rip- 
pling of the layers^^ and due to the disorder and electron- 
phonon coupling in the substrate 

In this paper we study the effects of randomness on 
the spin relaxation and spin injection. We show that 
spin relaxation reveals interesting quantum effects aris- 
ing from the non-commutativity of the momentum and 
coordinate-dependent randomness. The calculated spin 
injection can be observed in a wide range of frequencies, 
extended up to the electron Fermi energy. 

Model. We consider a two-dimensional electron gas 
with fluctuating Rashba SO interaction. In the absence 
of external fields, the Hamiltonian has the form (we use 
units with h = 1) H = Hq -\- Hso, where 
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H„^-— + Uir), 
Zm 



(1) 



Hso = -^a, {V„ A(r)} + '-ay {V„ A(r)} . (2) 

Here m is the electron effective mass, U (r) is the random 
potential leading to the momentum relaxation time Tp, 
and {,} stands for the anticommutator. The random 
Rashba field A(r) has zero expectation value (A(r)} = 
and correlation function 

(r - r') ^ (A(r) A(r')) - (A^) F (r - r') . (3) 
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where C\\{q) is the Fourier transform of the correlator 
Cxx (r) , (p is the angle between k and q, and q = k— k' is 
the momentum change due to spin flip scattering by fluc- 
tuations in SO field. Since the system is macroscopically 
isotropic in the xy plane, the coordinate-independent 
function Sk depends only on fc, that yields Sk = Sk', 
and thus we obtain 

StPi,= -%?, (9) 



FIG. 1: (Color online) Spin relaxation time as a function of 
the correlation radius R for different electron wave vectors k 
marked near the lines. 



The brackets (...) stands for the average over the disorder, 
and the range function F (r — r') depends on the disorder 
type. 

Spin relaxation. We begin with calculation of spin re- 
laxation rate due to random Hso interaction. The eigen- 
functions of Haniiltonian Hq (normalized to the unit 
area) are ~ e*'' '' Xa, where we included the spin in- 
dex a in the definition of momentum k = (k,cr), and Xa 
is the spin function. Matrix elements of SO interaction 

k|i?so|k ) are 
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{a \a^^ {ky + k'y) - ay (k^ + k'^)\ a') , (4) 



where A^-k' is the Fourier component of the random 
Rashba field. 

To calculate the spin relaxation time, we use the ki- 
netic equation for spin density matrix (see for instance 
Refs. [21,22,23]) 



(5) 



where rf is the spin relaxation time, 



Here, due to the absence of the regular contribution in 
the Hso, the commutator term in Eq.([5]) vanishes. There- 
fore, the entire effect of SO randomness is included in the 
collision integral^ 

St Pk - ^ E (2^id^'Pk'^k'k - ^'^'/'k - Pk%'^k'k) 

k' 

x5(ek - Ek') ,(6) 

with kinetic energy ek = k"^ /2m. We take in the 
form corresponding to the only nonzero z-spin compo- 
nent = -|- S'k , with being the equilibrium 
density matrix. The resulting macroscopic spin density: 
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Using Eqs. ^ and ^ we find 
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1 



m 
47r 



2k 



(10) 



(11) 



We employ the following form of C\\{q): 
CAA(q) = 27r(A2)i?2e-9fl 

where R is the length scale of variations in A. This form 
of correlator is realized when the Rashba SO coupling is 
formed by the z-component of electric field of random 
donors^^ symmetrically distributed on both sides of the 
QW at the distance L — R/2 from the QW symmetry 
plane. 

With Eqs. Uni) and dTT]) we obtain 
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x'^S^^'^ dx 
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[Ii{2kR) - Li{2kR)] 



(12) 



where Ii{x) and Li{x) are the Bessel and Struve func- 
tions, respectively, and 1/tso = m(A^). As a result, we 
obtain 



kR, 



kR-» 1, 



- = — M (13) 

'^k TT{kR)^/2, kR < 1. 

Equation ([13]) agrees with the results of Ref.[l3] for 
kR 3> 1, and shows that for given Tsq the relaxation rate 
rapidly decreases at small kR. Due to the anticommu- 
tator form of i?so, at small kR main contribution to the 
SO field comes from the derivatives of A(r). The increase 
in the relaxation time in this regime can be understood 
as a decrease in the disorder effect due to the averaging 
of _ffso over the area of 1/fc^. The dependence of the spin 
relaxation time on the correlation radius R is presented 
in Fig. 1. 

Combined resonance. Now we consider response of the 
system in a static magnetic field B to an external pe- 
riodic field A{t) = Aqc"*"*. The corresponding inter- 
action Vcxt = — (e/c)vA, where v = i[Ho + Hso,r] is 
the velocity operator, induces combined resonance caus- 
ing transfers between states with different spins and mo- 
menta. To study the spin dynamics, we retain only the 
spin-related part of the Hamiltonian and present it as: 

H = Ho + Hso+HB + V,^t, (14) 

A e 
Hb = -^{cr ■ n), 14xt = {(^xAy - ayA^) , (15) 




FIG. 2: (Color online) Possible spin-flip transitions. The ini- 
tial state is k' |, the final states are ki,2 i- Dashed circle 
corresponds to the transitions allowed by the energy conser- 
vation. The region in the momentum space with the size of 
the order of inverse correlation length 1/R, where the transi- 
tions can occur effectively, is also marked in the figure. If the 
final state is inside this area, the transition has a relatively 
high probability. 



where A — gfisB, g is the electron Lande factor, and n is 
the direction of B. We include the magnetic field via the 
Zeeman term, while neglect its orbital effects. The elec- 
tron energy spectrum is then spin-split, ekt.i = ek± A/2, 
where arrows correspond to the direction parallel and op- 
posite to the magnetic field. At realistic conditions, the 
splitting is much smaller than the chemical potential /.t 
of the degenerate electron gas. Here the periodic field 
leads to a disorder in V^xt due to the factor A(r). As a 
result, Vc-x_t causes transitions with the change in the elec- 
tron momentum and spin in a single process, as shown in 
Fig. 2. This is in contrast to conductivity, where the cou- 
pling of the external field to the disorder appears only 
through the disorder effect on the electron states, and 
the transitions are momentum-conserving. We calculate 
below the corresponding spin pumping rate. 

For given geometry of the external fields, denoted here 
as [g] the time evolution of the spin projected electron 

^(w, A) and spin- 



density riff' is due to the spin-gain J^^| 
loss (cj, A) processes: 

^^/l^l_(c.,A)-/il, 



.(^,A). 



(16) 



The concentration gain 
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[/(^k')-/M 5{s^-e^,-u.\ 
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is due to all possible transitions from occupied a' to 
unoccupied a states; a similar expression holds for the 
loss /o-^cr' (tj, A). The perturbation associated with the 
dipole moment acquired by electron spin in the presence 
of SO coupling^ has the form 
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(18) 




FIG. 3: (Color online) Spin pumping as a function of uj 
for different values of correlation radius R marked near 
the lines, (a) In-plane field and linearly polarized radia- 
tion. Weak peaks seen at small i? at = /i are mani- 
festations of the spin-split density of states. We use units 
{d{s)/dt)o = I/tt X (eVc^)^LMo. (b) 2-axis field with 
B and circularly polarized radiation. We use units 
(d(s) /dt)o = 2/7r X (eVc')Ag/r.o. 
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Substituting psp into p7p and averaging over spin dis- 
order, we obtain the component of the generation rate 



kk' 

X [/(^k') - fi'^k' + '^)] '^(^k - ^k' - ^) 



(20) 



where the coefficient K^^} is determined by the field 
configuration. Here we consider two geometries: (i) In- 
plane magnetic field B = (5,0,0), (sg' 



- {sx), lin- 
early polarized radiation Ao = {AoxjAoy) and (ii) z- 
axis magnetic field B — (0,0,5), {^Sg ^ — (s^), circu- 
larly polarized radiation Aq = Ao(x, ly). In the case 
(i) kI^}^^ = Aq,j.. The transition rates satisfy symme- 
try relation /cr'^cr('^, A) = /o-^cr' (w, — A), and dn\^^/dt 



[g] 







in Eq. (fT6|) vanishes if _B = 0. In the case (ii) 
for the transitions up-down and kI^}^^ — 2Aq for the 
opposite process. Therefore, for the circularly polarized 
radiation the spin pumping occurs even at i? — > 0. 

The exact formula for the pumping rate, valid in the 
general case /i > A/2, can be obtained after integrating 
Eg. ([20]) over the directions of q = k — k'. As a result one 
obtains (with y = kR) 



Due to charge conservation, d{n^a'^ + 7i'^})/dt — 0. Thus, 
pumping rate for the spin density component along the 
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^ J\y-y^,\ [47;2a;2 - (a;2 + 2m(a'A - u;)i?2)2]i/2 ' 

1/2 

where ymin,a-' = i? [max{0, 2m(/i — ct'A/2 — cj)}] , 
2/max,.' = R[2m{^i - ff'A/2)]i/2, y^, = [y2 ^ 2m(c^ - 
cr'A)i?2]i/2^ and I/tjo is the prefactor in Eq. (fH)) . 

For numerical calculations we use the following pa- 
rameters: electron effective mass for a Si/Ge (001) QW 
m = 0.19 mo (where mo is free electron mass), electron 
concentration per valley n — 5 x 10^^ cm, Fermi mo- 
mentum kp = 1.8 X 10^ cm^^, and the Fermi energy 
/i = 6.3 meV. To calculate d{sx) /dt we take magnetic 
field i? = 1 T leading to the spin splitting A = 0.12 meV 

for 5 = 2. The injection rates <i<^SB ^ /"^^ presented 
in Fig. 3. For both radiation polarizations, the peaks 
have the maximum position and the width on the or- 
der oi Lo fi/kpR. With the increase in R at given kp, 
the available momentum and energy ranges decrease, the 
peaks sharpen and shift to lower frequencies. For linearly 
polarized radiation (Fig. 3(a)), the pumping rate is linear 
in A. 



Conclusions We have studied spin relaxation and in- 
frared radiation-induced spin transitions in a 2D electron 
gas with the Rashba field disorder. Quantum effects re- 
lated to noncommutativity of the momentum and ran- 
dom Rashba potential lead to the decrease in the spin 
relaxation rate when the spatial scale of the randomness 
is smaller than the electron wavelength. In contrast to 
conductivity, external periodic electromagnetic field gen- 
erates a perturbation directly including the SO disorder, 
and, therefore, causing spin-flip accompanied by a mo- 
mentum change. As a result, electron spin density can 
be pumped by coupling of spins to the external peri- 
odic field in the frequency range up to the Fermi energy. 
These effects extend the abilities of manipulating spins 
in semiconductor structures. 

Acknowledgements. This work is partly supported by 
the FCT Grant PTDC/FIS/70843/2006 in Portugal and 
by Polish Ministry of Science and Higher Education as 
a research project in years 2007 ~ 2010. E. Sherman ac- 
knowledges support of the University of Basque Country 
UPV-EHU grant GIU07/40 and valuable discussion with 
M.M. Glazov. 



^ I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 

76, 323 (2004) 

^ E.I. Rashba and V.I. Sheka, Fizika Tverd. Tela (Leningrad) 
3, 2369 (1961) [Sov. Phys. Solid State 3, 1718 (1962)]; E.I. 
Rashba and V.I. Sheka, in: G. Landwehr and E.I. Rashba, 
Editors, Landau Level Spectroscopy, North-Holland, Ams- 
terdam (1991), p. 131 and references therein. 

^ E. I. Rashba and Al. L. Efros, Phys. Rev. Lett. 91, 126405 
(2003); Appl. Phys. Lett. 83 5295 (2003). 
M. Duckheim and D. Loss, Phys. Rev. Lett. 101, 226602 

(2008) . 

^ D. V. Khomitsky, Phys. Rev. B 77, 113313 (2008); Phys. 

Rev. B 79, 205401 (2009). 
^ K. C. Nowack, F. H. L. Koppens, Yu. V. Nazarov, L. M. 

K. Vandersypen, Science 318, 1430 (2007). 

M. Pioro-Ladriere, T. Obata, Y. Tokura, Y.-S. Shin, T. 

Kubo, K. Yoshida, T. Taniyama, S. Tarucha, Nature 

Physics 4, 776 (2008). 
« Y. J. Bao and S. Q. Shen, Phys. Rev. B 76, 045313 (2007). 
^ E. Ya. Sherman, Phys. Rev. B 67, 161303(R) (2003). 
^° M. M. Glazov and E. Ya. Sherman, Phys. Rev. B 71, 

241312(R) (2005). If the Rashba field is produced by the 

random donors, (A^) behaves as and, therefore, the 

spin relaxation time decreases with the increase in R. 

L. E. Golub and E. L. Ivchenko, Phys. Rev. B bf 69, 

115333 (2004) 

C. P. Moca, D. C. Marinescu, and S. Filip, Phys. Rev. B 

77, 193302 (2008). 

Y. Tserkovnyak and S. Akhanjee, Phys. Rev. B 79, 085114 

(2009) . 



Z. Wilamowski, W. Jantsch, H. Malissa, and U. Rossler, 
Phys. Rev. B 66, 195315 (2002). 
^5 M.-H. Liu, K.-W. Chen, S.-H. Chen, and C.-R. Chang 
Phys. Rev. B 74, 235322 (2006). 

J. D. Koraiek, C. P. Weber, J. Orenstein, B. A. Bernevig, 
Shou-Cheng Zhang, S. Mack, and D. D. Awschalom, Na- 
ture 458, 610 (2009). 

G. M. Miiller, M.l Romer, D. Schuh, W. Wegscheider, J. 
Hiibner, and M. Oestreich, Phys. Rev. Lett. 101, 206601 

(2008) . 

i« V. V. Bel'kov, P. Olbrich, S. A. Tarasenko, D. Schuh, W. 

Wegscheider, T. Korn, C. Schiiller, D. Weiss, W. Prettl, 

and S. D. Ganichev, Phys. Rev. Lett. 100, 176806 (2008). 

D. Huertas-Hernando, F. Guinea, and A. Brataas, Eur. 

Phys. J. Special Topics 148, 177 (2007). 
^° C. Ertler, S. Konschuh, M. Gmitra, and J. Fabian, Phys. 

Rev. B 80, 041405 (2009). 
2^ S.A. Tarasenko, Pis'ma v ZhETF 84, 233 (2006) [JETP 

Letters 84, 199 (2006)]. 
^2 D. Culcer and R. Winkler, Phys. Rev. B, 76, 195204 

(2007) ; D. Culcer, preprint cond-mat arXiv:0904.1999] 

(2009) ; R. Winkler, D. Culcer, S. J. Papadakis, B. Habib, 
and M Shayegan, Semicond. Sci. Technol. 23, 114017 

(2008) . 

M. M. Glazov, Sol. State Commun. 142, 531 (2007); Phys. 
Rev. B 70, 195314 (2004). 

T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 
437 (1982). 



